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^ ■ Abstract 

The discrepancies in claims of experimental evidence in the search for weakly in- 
teracting massive particle (WIMP) dark matter necessitate a model for ionization 
_ efficiency (the quenching factor) at energies below 10 keV. We have carefully studied 

O ; the physics processes that contribute to the ionization efficiency through stopping 

' . power. The focus of this work is the construction of a model for the ionization ef- 

\-^ \ ficiency in germanium by analyzing the components of stopping power, specifically 

c/3 ' that of the nuclear stopping power, at low energies. We find a fraction of the ZBL 

nuclear stopping power can contribute to ionization efficiency. We propose a model 
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that corrects the missing contribution to ionization efficiency from the ZBL nuclear 
stopping power. The proposed model is compared to previous measurements of ion- 



C^ ■ ization efficiency in germanium as well as that of other theoretical models. Using 

J^ , this new model, the thresholds of both CDMS II and CoGeNT are analyzed and 

.^ I compared in terms of the nuclear recoil energy. 
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1 Introduction 



In the search for dark matter through direct detection of weakly interact- 
ing massive particles (WIMPs), it is important to create well defined models 
for nuclear interactions at low energies. WIMPs are believed to only interact 
via weak force and gravity; hence the interaction cross section with ordinary 
matter is extremely small. If they are to be detected, it will be from their 
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scatterings off nuclei in the target material. Measuring the recoil of the nu- 
clei (directly or indirectly) will give an estimate of the mass of the incoming 
WIMPs. 

DAMA/Nal [I] and DAM A/LIBRA [2] have observed a model independent, 
annually modulated signal that has been interpreted as evidence for dark mat- 
ter. This claim has been met with skepticism within the global dark matter 
community. Before the end of 2010, CDMS [3], EDELWEISS g], XENONIO [5], 
XenonlOO [B], and CoGeNT [7] published limits that are clearly incompati- 
ble with the DAMA signal region [TP|5] by assuming coherent scalar WIMP 
interactions. Moreover, the combination of the upper limit contours from the 
CDMS II [9], CRESST I [10], PICASSO [IT], NAIAD [12], ZEPLIN I [13], 
EDELWEISS [H], SIMPLE [15] and Super-Kamiokande [E] experiments in- 
dicates an inconsistency when interpreting the DAMA signal in terms of a 
WIMP-nucleon spin-dependent coupling within the standard halo model. A 
global analysis by Fairbairn and Schwetz [17] and Savage et al. [18] disfavors 
the DAMA results as a WIMP-nucleon interaction. 



However, the recent CoGeNT annual modulation signal [IH], released in June 
2011, is a star of excitement for WIMPs in a low mass range of 7 GeV. In 
response to this excitement, CDMS II [20] performed a low threshold analysis 
and found there is no signal in the region of interest that is comparable with 
CoGeNT. On March 5th of 2012, CDMS II released the analysis of annual 
modulation in the low energy CDMS II data [21] and found no evidence of 
annual modulation. Note that both CoGeNT and CDMS II are germanium- 
type experiments. This discrepancy signifies the importance of understanding 
the detector response to low energy nuclear recoils. 

Since CDMS II has n-7 discrimination capability down to the low energy 
region, and CoGeNT and DAMA do not posses such ability, the possibil- 
ity that dark matter particles interact predominantly with electrons as op- 
posed to nuclei, discussed by many authors including the DAMA collabora- 
tion [22|23f24ll25f26] . is an interesting approach that is under consideration. 
This requires some extension of the Standard Model but nothing forbids such 
a possibility. However, this possibility is now fading because of the most recent 
analysis of electronic recoil events in CDMS II low energy data j^ with no 
evidence of annual modulation. In addition, the result from CRESST-II |27j . 
which claimed the excess of events in their energy range that indicates a 
WIMP mass of 13 GeV, differs from CoGeNT's mass of WIMPs. CRESST- 
II does possess a good n-7 discrimination in their energy region of interest. 
This is another indication of the importance of understanding the detector re- 
sponse to low energy nuclear recoils. Note that XenonlOO shows no evidence of 
WIMP mass in the region of 7 GeV to 13 GeV [^. There is also an alternative 
approach to the explanation of experimental results assuming an 0-helium 
universe [251. 



Lack of a coherent picture to explain the results from DAMA, CoGeNT, and 
CRESST-II requires more attention be paid to the interpretation of the exper- 
imental results, such as the ionization efficiency, which is critical to the detec- 
tion threshold of energy, hence the region of interest and the WIMP masses. 
Since CoGeNT and CDMS II are both germanium detectors and the results 
conflict with each other, it is necessary to study the ionization efficiency of 
germanium detectors. Several theoretical models exist for ionization efficiency, 
however the best model for interpreting the experimental data has yet to be 
established. Many experiments adopt a model proposed by Lindhard et al. in 
the 1960s [21] with different values of a constant, fc, in their data analysis. 
For example, within the limits of A; = 0.1 io k = 0.2, k = 0.15 [30]|31II32||33] . 
k = 0.157 [31135] . k = 0.159 [36] and A; = 0.2 [301I37J were used in the inter- 
pretation of the experimental data. The difference in the ionization efficiency 
given by different values of k can be as large as 30%. Without knowing which 
value of k is correct in terms of physics in the low energy region, it is diffi- 
cult to justify the threshold of nuclear recoil energy and the mass of WIMPs. 
Therefore, a model that demonstrates physics processes in the low energy re- 
gion is needed. This model will also justify which value of k should be used in 
Lindhard's model. The model presented in this paper has been developed by 
re-examining the components of stopping power, especially that of the nuclear 
stopping power, at low energies. 



2 Developing the Model for Ionization Efficiency 



Previous work by Mei et. al. [38] developed models for the quenching factor in 
noble liquids. This work modifies the previous work to look at the ionization 
efficiency in germanium, taking another look at the stopping power of the 
germanium ions produced in the WIMP-nucleus interaction in germanium. 



2.1 Electronic Stopping Power 



The electronic stopping power is well understood as the orbital electrons of the 
surrounding atoms that are ionized by inelastic collisions with the recoiling ion 
along its path. This portion of energy loss by a recoiling nucleus contributes 
to the ionization efficiency at a level of 100%. To calculate the electronic stop- 
ping power, the heavy ion scaling rule was applied to calculate the stopping 
power of germanium ions in germanium from the stopping power of protons 
in germanium at the same velocity: 

2, 



Sce — {(.Zce) Sp. 



The effective charge fraction, ^, was approximated by Brand and Kitagawa [39] 
and modified by ZBL [IQ] as 



^ = g + 0.5-(l-g)— -/n 



4A-2.02 
1.919 



(2) 



where vp is the Fermi velocity, vq the Bohr velocity, A is given by 



^ 0.48(1 -g)3 
A = Y ao 
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(3) 



with q given by 

q = l-exp f0.03732922/°-3 + 0.0571791y°-^ - 1.09793?/^ + 0.072649?/2l , (4) 



and ao is the Bohr radius, i/r is a dimensionless quantity related to the relative 
velocity by 



Vr 



^0 • ^Ge 



(5) 



The relative velocity is given in [41], and is related to vp and the ion velocity, 
wi, by 



l^F(l + l^^-^,,^J,iiv^<vp, 



^1 1 + 



5^ 



if vi > vp- 



(6) 



Figure [T] shows the calculated electronic stopping power in germanium. 



2.2 Nuclear Stopping Power 



Nuclear stopping power is less understood at the low energies induced by the 
recoiling ion in a WIMP collision. Nuclear stopping represents elastic collisions 
between the recoiling ion and the surrounding atoms in the target. If one knows 
the form of the repulsive potential between them, it is possible to calculate 
the nuclear stopping power. The general calculation of nuclear stopping power 
is performed with the following equations from ZBL 
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Fig. 1. Shown is the calculated electronic stopping power in germanium for the low 
energy region of nuclear recoils. 

where Mi and M2 are the projectile and target masses (amu), and Zi and Z2 
are the projectile and target atomic numbers, respectively, e is the reduced 
energy given by 



32.53M2E 



ZiZ2(Mi + M2)(Z0-23 + z0- 
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and the reduced nuclear stopping power 




;n.(l+1.1383£) 



2(€+0.01321eO-2i226_|_o, 19593^0. 5) 

;n(£) 
2e ' 



oTsT, fore < 30, 



fore > 30. 



(9) 



Note E is kinetic energy of projectile in keV. 

We show the calculated nuclear stopping power in Figure [2] using the above 
formulas. 



In general, this part of energy loss does not contribute to ionization efficiency. 
However, the standard deviation of the fit of the universal ZBL repulsive 
potential to the theoretically calculated potentials is 18% above 2 eV [ID] . 
This indicates that nuclear stopping power calculated with the ZBL repul- 
sive potential possesses large uncertainty. Comparing with a more accurate 
repulsive potential generated from self-consistent total energy calculations, 
the calculation of nuclear stopping power using ZBL repulsive potential can- 
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Fig. 2. Shown is the calculated nuclear stopping power in germanium for the low 
energy region of nuclear recoils. 



not be treated totally as the nonionizing nuclear stopping power. Consider- 
ing density-functional theory, the local-density approximation (LDA) for elec- 
tronic exchange and correlating both the electronic part and the inter-nuclear 
Coulomb repulsion as a function of the distance between atoms in a dimer 
bond, it became obvious to us that a correction needed to be derived to ac- 
count for the fraction that contributes to ionization efficiency |12P3] . 



2.3 Correction Factor 



When atoms receive significant recoil energies, they will be removed from their 
lattice positions, and produce a cascade of further collisions in the lattice. In 
the low energy range, the most prominent component of nuclear stopping 
power is the non-ionizing nuclear stopping power that describes the rate of 
energy loss due to atomic displacements as a recoil nucleus traverses a ma- 
terial |^45|46] . In the calculated nuclear stopping power using the ZBL re- 
pulsive potential, there is a small fraction of the nuclear stopping power that 
contributes to the calculation of ionization efficiency. It is the goal of this 
model to incorporate this fraction to develop a more accurate model of the 
ionization efficiency. The calculation of the fraction of the ZBL nuclear stop- 
ping power that contributes to the ionization efficiency is performed by the 
calculation of the non-ionizing energy-loss (NIEL) using the Wentzel-Moliere 
differential cross section discussed in Section 2 [Eq. (15)] of reference [Hj. The 



NIEL can be calculated with 
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where E is the kinetic energy of the incoming particle (in keV), E^ is recoil en- 
ergy (the kinetic energy transferred to the target atom) of the atom (in keV), 
and L(i?r) is a fraction of the recoil energy, E^, which undergoes displacement 
processes. The expression of L(E'r) is often called Lindhard partition function, 
which can be found in Ref. [771^8] and in Equations (4.94, 4.96) of Section 
4.2.1.1 in Ref . ^9] . Er{de) = ErL{Er) is the damage energy, i.e., the energy 
deposited by a recoil nucleus with kinetic energy Er undergoing displacement 
damage inside the medium. The integration lower limit of Er{d) = 23 eV in 
germanium [50] is the energy necessary to displace the atom from its lattice 
position. The upper limit of Ej.{max) is the maximum energy that can be 
transferred during a single collision process. Taking the difference between 
Eq. Hand Eq. [10] divided by the nuclear stopping power (Eq. [7]) gives the frac- 
tion of the nuclear stopping power that contributes to the ionization efficiency. 
This fraction of the ZBL nuclear stopping power was found to be 

Cf = 6.2 X 10-2E°■l^ (11) 



where Er is the recoil energy of the target atom received from the incoming 
particle. Figure [3] shows the fraction of the ZBL nuclear stopping power that 
contributes to the ionization efficiency. 



2.4 Ionization Efficiency 



After calculating the fraction of the ZBL nuclear stopping power that con- 
tributes to the ionization efficiency, it is applied to the ionization efficiency in 
the following way: 



where r]{Er) is the electronic stopping power, p{Ej) is the ZBL nuclear stop- 
ping power and Cf is the correction fraction. 

With this fraction of the ZBL nuclear stopping power included in the ioniza- 
tion efficiency, the theoretical model will be a more accurate picture of the 
experimental signal. Figure H] shows the calculated ionization efficiency as a 
function of nuclear recoil energy. 
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Fig. 3. The fraction of the ZBL nuclear stopping power contributing to the ioniza- 
tion efficiency as a function of nuclear recoil energy. 
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Fig. 4. The calculated ionization efficiency as a function of nuclear recoil energy. 
The fitted function can be expressed as: 
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This equation is valid for recoil energy, E^, greater than 1 keV and less than 
100 keV. The corresponding visible energy, Ey, can be obtained by multiplying 



the ionization efficiency by the recoil energy. 

It was important to ascertain whether nuclear recoils off the naturally occur- 
ring isotopes of germanium in the detectors could produce large differences 
in the ionization efficiency. Calculations were made for each natural isotope 
of germanium (70, 72, 73, 74, 76) and it was determined that the spread in 
ionization efficiency was less than 2%. This will contribute a difference of no 
more than a ±0.04 keV change to the visible energy for a recoil energy of 2 
keV. Therefore, the average atomic mass was used in all calculations. 



3 Comparison to the Experimental Data 



Previous works have measured the relation between ionization efficiency or vis- 
ible energy and the recoil energy in germanium [5I]]|31|32|33|34|36|37|51|52|53] . 
The calculated ionization efficiency in this work is compared to the measured 
data points from various measurements in Figure [5l As can be seen in Fig- 
ure [5l the measured data points from different experiments do not entirely 
agree with each other. This indicates the systematic errors involved in differ- 
ent measurements. 

Depending on the experimental techniques, the systematic errors can be dom- 
inated by different sources. For instance, thermal neutrons can be used to 
measure nuclear recoils down to sub-keV region. In this measurement, unless 
the out-going gamma rays are fully absorbed and measured by another de- 
tector in coincidence, Compton scattering of the out-going gamma-ray inside 
the germanium detector can contaminate the visible energy. In the case of 
measuring nuclear recoils with elastic scattering, the multiple scatters must 
be excluded by a Monte Carlo simulation. In addition, the scattering angle 
and the time of flight of the out-going neutrons have to be measured very 
accurately. However, there is a limitation on how fast the electronics can pos- 
sibly respond to the signal. If one uses inelastic scattering to measure nuclear 
recoils, the de-excitation of gamma rays through Compton scattering in the 
detector can contaminate the signal. Therefore, all of the measurements must 
be implemented together with a reliable Monte Carlo simulation, which is well 
calibrated by gamma rays in the same energy range, to reduce the systematic 
error and deliver reliable ionization efficiency. This is to say that a best fit to 
some of data points in the energy region of interest will not necessarily rep- 
resent the true ionization efficiency. The physics that involves the ionization 
efficiency must be understood. 

Because of a large spreading of experimental data points, the predicted ion- 
ization efficiency from the model agrees with some data points but disagrees 
with other data points, with the majority of measurements above the model. 
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Fig. 5. A comparison between the calculated and the measured ionization efficien- 
cies. Note that the measured data points are from Jones et al. |30|51j . Barbeau et. 
al. [37], TEXONO [52], CDMS [20], Messous et al. [M], Chasman et al. [3l], Shutt 
et al. [32], Sattler et al. [36], Baudis et al. [33], and Simon et al. [53] . 

The CDMS II data was collected comparing the ionization yield to the recoil 
energy [20]. The possible discrepancy could be due to the lack of the calibration 
of the phonons in the low energy region. 



4 Comparison to the Previous Models 



The most accepted model for calculating the ionization efficiency, e, 
element is developed by Lindhard et. al. [29] 



l + k-g{e) 



of any 



(14) 



where e = 11.5ErZ ^/^ for a given atomic number, Z, k = 0.133Z'^'^A ^/^, 
and c/(e) is fitted by ^(e) = 3e°-i5 + 0.7e°-6 + e. 

Often two different values of the scaling constant have been used in conjunction 
with the experimental data, one that agrees better at the higher recoil energies 
and one at the lower. These are plotted against the model for k = 0.1 and 
k = 0.2. The calculated value of k for germanium is k = 0.159. This is plotted 
with the traditional k values in Figure [61 It is interesting to note that at low 
energies k = 0.159 agrees well with the model proposed by this work. Thus, 



10 



if Lindhard's model is to be used, the proper value of k for low energies is, 
k = 0.159, for germanium. Also note that at higher energies the proposed 
model tends to agree with the Lindhard model with k = 0.1. 

A best fit to the data between 1 keV and 10 keV has been used by J.I. Collar 
et al. PTIIMJ and others [21133] in order to interpret the CoGeNT and CDMS 
II results: 



E„ = 0.199K 



1.12 



(15) 



This best fit appears to follow the Lindhard model k = 0.2 at low energies 
and deviate to become parallel with Lindhard A; = 0.1 as the energy increases. 

Figure [6] shows a comparison between the proposed model, Lindhard's model, 
and Collar's best fit. 
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Fig. 6. Visible energy versus nuclear recoil energy from different models. 



5 Validation of Models with Experimental Data 



Both CDMS II and CoGeNT are low threshold detectors accessible to low en- 
ergy recoils induced by low mass WIMPs. Therefore, the ionization efficiency 
in the low energy region is critical to both experiments. To validate a reliable 
model for ionization efficiency in the low energy region, we plot the experimen- 
tal data points together with several theoretical models including the proposed 
model in this work in Figure [71 It is quite clear in Figure [7] that Lindhard's 
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Fig. 7. Ionization efficiency as a function of energy of nuclear recoils. Shown is a 
comparison between the measured data points and the popular models including 
the proposed model in this work. 

model with k = 0.1 ot k = 0.2 does not agree with the measured data points in 
particular for the low energy region. This is because k = 0.159 for germanium 
according to Lindhard's theory [29]. There is no theoretical ground to choose 
k = 0.1 or k = 0.2 in Lindhard's model. li k = 0.159 is chosen for Lindhard's 
model, there is a fair agreement between Lindhard's model and the proposed 
model in this paper. The difference is on the average of 10% across the energy 
range from 1 to 100 keV. This can be understood according to Lindhard's the- 
ory in which Lindhard pointed out that his model possesses uncertainty in the 
low energy region [56]. It is thus important to notice that k = 0.159 must be 
chosen if one uses Lindhard's model for germanium detector. The 10% error 
could shift the analysis threshold by ±0.2 keV at the recoil energy of 2 keV. 
This is a large enough error to generate argument and result in discrepancy 
in the analysis. A numerical comparison between different values of constants 
for Lindhard's model, the best fit of Collar et ai, and the proposed model is 
shown in Tabled] As can be seen in Figure [7] and Tabled! the proposed model 
is reliable from 1 keV to 100 keV. 



Application of the Proposed Model to the Thresholds of CDMS 
II and CoGeNT 



CoGeNT claimed the evidence of annual modulation indicating a mass of ~7 
GeV WIMPs dSJ. CDMS II failed to confirm CoGeNT's claim [ZH. Both ex- 
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Table 1 

Comparison between visible energies at a given recoil energy for the models and 
the Collar's best fit (all energies in keV). 









Lindhard 


Recoil Energy 


Proposed Model 


Collar 


k = 0.1 


k = 0.159 
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0.169 


0.199 
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0.361 
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periments use germanium detectors located in the Soudan Mine. CoGeNT 
cannot discriminate electronic recoils from nuclear recoils while CDMS II pos- 
sesses a good capability of discriminating n-7 events. The discrepancy between 
two experiments lies in whether CoGeNT sees nuclear recoil events that are 
below the detection threshold of GDMS II. 



CDMS II used a 2 keV nuclear recoil threshold in the analysis of the low 
energy region [20] in which the ionization efficiency is 15.8% [20]. We calculate 
a difference in the energy scale by °'^^^°^~^V^^ = 12.5%, where 18.05% is the 



0.1805 



ionization efficiency at 2 keV nuclear recoil energy from the proposed model. 
To interpret the CDMS II threshold with the proposed model, we multiply it 
by the ratio of the ionization efficiencies: 2 keVx^4ll = 1-75 keV. This is the 



0.185 



CDMS II threshold in the energy scale of the proposed model assuming that 
the CDMS II energy scale determined by the phonon signal is correct. 



Using the proposed model, the analysis threshold of CoGeNT, 0.5 keV elec- 
tronic equivalent energy, is calculated to be 2.7 keV nuclear recoil energy. The 
lowest threshold of CoGeNT, 0.4 keV electronic equivalent energy, corresponds 
to 2.23 keV of nuclear recoil energy. Thus, the CDMS II threshold of 1.75 keV 
is lower than both the 2.7 keV analysis threshold of CoGeNT and the 2.23 
keV lower limit of the CoGeNT detector. Figure M shows a threshold energy 
comparison between CoGeNT and CDMS II. As can be seen, the CoGeNT is 
fully contained by CDMS II using our model in this paper. 

Using Collar's best fit equation, Eq. [151 the lowest electronic equivalent energy 
of 0.4 keV corresponds to 1.87 keV of nuclear recoil. The difference between 
the Collar's best fit equation and the model proposed in this work is about 
18% in calculating nuclear recoil energy for a visible energy of 0.4 keV. 
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Fig. 8. The threshold energy of CoGeNT versus that of CDMS II for WIMPs with 
mass of 7 GeV and an interaction cross section of 2.5 xlO~^^ cm^. Included in this 
plot is the event rate verses recoil energy of other common (or considered) targets 
for WIMP dark matter detection. 

7 Conclusion 



We have developed a new model for ionization efficiency by analyzing the 
components of stopping power and the fundamental physics that must be 
considered to accurately understand low energy nuclear stopping power. This 
model is ffist calculated for germanium as it is the simplest of detectors. The 
result is compared to both experimental data and previous theoretical models. 
When compared with the experimental data points, this model passes in the 
middle of the data points. It shows some agreement with a few sets of data 
but discrepancies between others. Many sets of experimental data did not 
span the entire range needed, especially at low energies data points are scarce, 
and so verification will be necessary as a next step. When compared to the 
theoretical models, there is a fair agreement between Lindhard model with 
k = 0.159 and our model at low energies. This value of k is the calculated 
value for germanium [Z = 32, A = 72.64) in the Lindhard theory. 

The analysis threshold of the CoGeNT detector is 2.7 keV nuclear recoil energy 
using the proposed model in this paper, which is higher than that of the 1.75 
keV threshold of the CDMS II detector in the same model. Thus CoGeNTs 
region of interest in completely contained within that of the CDMS II detector. 
Based on the analytic model developed here and the published results of the 
two collaborations [T9|20] . the excess events above the background seen by 
CoGeNT seem to not be nuclear recoils. 
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Work will continue on the proposed model to verify its accuracy in the range of 
1 keV to 10 keV. Understanding the ionization efficiency in germanium is the 
first step to building a set of models for ionization efficiency in other materials 
employed by dark matter detectors. This will provide reliable analytical models 
for the interpretation of experimentally observed signals. 
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